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[3] Modeled by a strong wall repulsion. 

[4] Positivity rests on numerical evidence [14]. 

[5] This is usually also assumed in dealing with the ergodic hypothesis, or more generally 
with principles which cannot be assumed to be independent from the basic laws of 
motion. 

[6] The assumptions (C) may just be too strong and/or difficult to verify for the model in 
Eq.(l). Their strength can be seen from the fact that they imply (R) above, [29], [26], 
and (hence) the ergodic hypothesis if 7 = 0, [29], [26]. The density in C could be replaced 
by the requirement of density in A alone; this would be more general without affecting 
our conclusions. Furthermore our main conclusions could still be reached under far 
weaker assumptions, as we think that the consequences of (C) relevant for our analysis 
naturally follow in the frame of the theory of singular hyperbolic systems of [24] . 

[7] We identify this with the actual time interval between them, neglecting its fluctuations 
(which could be easily taken into account leading to the same end result). 

[8] This is just as in the Ising model where one cannot compute correctly the thermody- 
namic limit average of the magnetization (or the average of any extensive quantity) for 
a finite volume r by using the finite volume Gibbs distribution with the same volume r 
(without incurring large, size dependent, errors). One must use a larger volume 6 ^> r, 
except in one dimension where it is well known [26] that the probability distribution of 
the magnetization would be off "only" by a factor bounded above and below for each 
spin configuration by a r-independent constant. Our case is very similar, as follows 
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to the angles between the manifolds at S'~ T / 2 x and at S T ^ 2 x. But assumption (C) implies 
transversality of their intersections, so that the ratio between JF t (x) and i*V(x) is bounded 
away from and +00 by (x, r)-independent constants. Therefore the ratio Eq.(4) is equal 
to e p( 2JV - 2 X <7 ) t o^ as we n as Iq that found in [15], up to a factor bounded above and below by 
a (r,p) independent constant, i.e. to leading order as r — > 00. 

We note that this result can be considered as a large deviation result (both in N and 
tor): its peculiarity is the linearity in p. On general grounds one might expect that the 
deviations of ( 2 jv-2)t T ^°S"^"( X )' or °f Pi fr° m its average value 1 have a probability density 
Tr(p) expressible in terms of some "free energy" function ((p) as tt (p) = e -C(p)rto_ Here ((p) 
describes the multifractal nature of the observable T T : a "non multifractal" distribution 
would correspond to a delta function distribution e^ p ^ Tt ° = 6(p — 1). Noting that e(p) = 
C(p) ~ C(~p) = ~^r\°E 7r (p)/ 7r (~P) i s an °dd function of p, we expect that e(p) = c (a) (p + 
■S3P 3 + s 5 p 5 + . . .) with c > and Sj ^ 0, since there is no reason, a priori, to expect a 
"simple" (i.e. with linear odd part) multifractal [9]. Thus p-linearity (i.e. sj = 0) is a key 
test of the theory and a quite unexpected result from the latter viewpoint. Recall, however, 
that as noted in remark (ii) after Eq.(3) the exponent (2N — 2)t T(cr)p is correct up to 
terms of 0(1) in r (i.e. deviations at small p must be expected). This shows that Ruelle's 
principle can indeed be tested on many particle systems in statistical mechanics and in fact 
is in agreement with the computer results obtained in [15]. The same theory would also 
apply to other models, like the ones in [11] or in the time reversible case in [12] for which, 
to date, no corresponding experimental results exist, yet. 

A similar test might even apply to fluid dynamics, for which (R) was originally devised. 
In fact the above linear multifractality (of the odd part of ((p)) might be observable in 
the high turbulence regime in fluid mechanics at least in numerical experiments for models 
in which no friction acts in the inertial range. The fluid is then described in the inertial 
range by the reversible Euler equations, while the dissipation below the Kolmogorov scale 
can be modeled by a gaussian thermostat thus making the equations describing the system 



8 



The main experimental result corresponds in the present context to the computation of the 
logarithm of the ratio of the probability that a T (x) = p to that of a T (x) = —p. The result, 
fig. 2 of [15], is that this quantity as a function of p is a precise straight line with slope 
2NTt (<j) for r large. 

We now use the very formal expression of p defined via Eq.(3) to study some statistical 
properties of JF t (x) and compare the result with [15]. We argue that the ratio of the 
probability of a T (x) £ [p } p-\- dp] to that of a T (x) £ [— p, — p-\- dp] is, using the notations and 
the approximation to p in Eq.(3): 

Sj, q r (xj)=p ^-m,t( X j) ^ 
Sj, a r (xj) = -p ^-«,r( X i) 

We evaluate Eq.(4) by establishing a one to one correspondence between addends in the 
numerator and in the denominator, aiming at showing that corresponding (i.e. with the 
same j) addends have a constant ratio which will, therefore, be the value of the ratio in 
Eq.(4). This is made possible by the time reversal symmetry which is the extra information 
we have with respect to [29], [10], [26]. In fact the time reversal symmetry (B) can be shown 
to imply that with £ also i£ is a MP. Since the intersections of MP are still MP we can assume 
that £, hence also will be time reversal symmetric, i.e. that for each j there is a j' such 
that iEj = Eji. By using the identities S~ T (S T x.) = x, and S~ r (iS~ T x.) = ix. (time reversal) 
and iW™ = Wf x , one can deduce, with A SjT (x) defined after Eq.(2): a T (x) = — a T (zx) and 
A„ jT (zx) = A S)T (x) _1 . The ratio Eq.(4) can therefore be rewritten as: 

Sgj,a r (xj)=p ^-m,t( X j) _ Sgj,a r (xj)=p ^m,t( X j) ^ 
S£jj,a r (xj) = -P ^-«,r( X j) S£Jj,a r (xj)=p A s , T (Xj) 

Then the ratios between corresponding terms in Eq.(5) are equal to F T (x.j) = A~^.(xj) • 
Ajl(xj). This is almost ^v(xj) = e a ^( x j) ( 2JV - 2 )( <7 )*o T which is j -independent (because 
a T (xj) = p). In fact, the latter is the reciprocal of the determinant of the jacobian ma- 
trix of S } i.e. the reciprocal of the total phase space volume variation, while F T (x.j) is only 
the reciprocal of the product of the variations of two surface elements tangent to the sta- 
ble and to the unstable manifold in Xj. Hence i*V(x) and JF t (x) differ by a factor related 
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/ u Tj t(^x)G'(x) 



(3) 



3) Consider the limit as r, T — > oo, with T — > oo fast enough compared to r so that in each 



parallelogram the weights in Eq.(3) have a small relative variation. 

It follows then from Sinai's work [29] and [10], [26] that the limit exists and is the 
statistics p of the Liouville distribution [21] and that p does indeed verify (R). The Ej may 
be interpreted as the dynamical states of [15]. The p T ^ in Eq.(3) plays the role of the finite 
volume Gibbs distributions in statistical mechanics. 

We now turn to the numerical experiment in [15]: 
(i) the main results concern properties of the entropy production over a generic time interval 
rt , during which the trajectory moves between S'~ T / 2 x and S T ^ 2 x, if t is the average time 
interval between two successive collisions, [7] (x is the middle point of this trajectory seg- 
ment). The entropy production is defined here, following [23], by ~ 2)a(S J x.). 
It will be convenient to rewrite the latter expression by setting (2iV — 2)rto(cr) T (x) = 
(2iV — 2)rt (cr) a T (x), where (<t) t denotes the time average of the entropy production rate a 
over the time interval rt . This defines a fluctuation variable a T (x) with (forward) average 
equal to 1 (because the infinite time average of (<r) T is the above defined (cr), see (A)). The 
exponential JF t (x) of the entropy production is the reciprocal of the phase space contraction 



(ii) For such an observable, which is strongly r dependent, it may be doubtful to use p T ^ to 
estimate the probabilities relative to //, even if r, T are very large. [8] However, it is possible 
to prove, if (C) holds, ( [21]), that not only the average of T T but even the probabilities of the 
various values of T T can be computed by using p T ^ of Eq.(3). The error on the probability 
attributed to each individual Ej would then consist of a factor bounded, above and below, 
by t -independent positive constant. 

In the computer experiment [15] one measures the entropy production log jF T (x), as seen 
on a stretch of time r short compared to the experiment duration T, repeatedly T/r times 
(see also [18]). This amounts to studying the p, distribution of the fluctuation a T (x), in (i). 



itself. 
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Most results deriving from (C) really come from the fact that (C) implies the existence 
of Markov partitions (MP), which permit us to give a rather simple description of the dis- 
tribution Ji, (the SRB, see Eq.(3) below), upon which all our deductions are based. The 
MP are tilings £ of C with suitable small sets E called "parallelograms", naturally related 
to the stable and unstable manifolds (with their "axes" parallel to the stable and to the 
unstable manifolds), see [10], [26] for a complete description. A property that always holds 
is the time invariance of the MP: if £ is a MP also the partition S£ obtained by evolving 
in time the elements of £ is such. A key property of MP is that the partition obtained by 
considering the intersections of all the tiles of two MP is still a MP. Hence, given a MP 
£, one can construct other much finer ones: a typical method is to consider the partitions 
S q £ and to define the partition £j obtained by intersecting all the tiles of all the partitions 
S q £ with —T<q<T. The parallelograms of the finer partitions may become as small 
as desired by taking T large (by the hyperbolicity). The description of the statistics p, of 
the Liouville distribution, in terms of the family of finer partitions associated with £, is 
done as follows: 

1) Given two times r and T, let Ej be the parallelograms of £j labeled by j, and let Xj be a 
point in Ej (arbitrarily fixed); define a probability distribution by attributing to each 
Ej £ £t a weight: 

Kfa) = U K'iS^) (2) 

k=-r/2 

where A„(x) is the absolute value of the determinant of the jacobian matrix of S as a map 
of to Wg x , so that the weight A~ 1 (x J ) is the inverse of the expansion coefficient of the 
map S T evaluated at S~ T ^ 2 Xj, i.e. at the initial point of a motion spending half of the 
time "before" reaching Xj and half "after" . For later use we define corresponding quantities 
associated with the stable manifold, denoted by A s (x) and A SjT (x). 

2) After normalization the above weigths define a distribution by requiring that the 
average of any (smooth) function G with respect to the distribution will be: 
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t — > x(t) is a solution of (1), so is z(x(— t)); [16,15]. 

We note that (A) implies the existence of an invariant set A, which we will call the attrac- 
tor } of zero Liouville probability (if 7 ^ 0) but probability 1 with respect to the asymptotic 
statistics p of the motions with random initial data (with the Liouville distribution), so 
that p(A) = 1. Property (B) implies a direct relation between the statistical properties of 
the forward it — > 00) and backward it — > —00) motions, although described in general by 
different statistics: e.g. they have the same set of Lyapunov exponents. 

To apply Ruelle's principle, in the spirit discussed above, we assume that, in view of 
the particle collisions, the behavior on A is as if A were hyperbolic. There may well be 
corrections to this: proceeding by ignoring such possibilities is our interpretation of (R). In 
practice this means that we suppose that the corrections become negligible as N — > 00. [5] 
The strongest form of the hyperbolicity assumption is: 

(C) (a) at every point x of C one can define stable and unstable manifolds W£ } W£ dense 
in C, transversal and covariant (i.e. they form an angle bounded away from 0, it when they 
cross and SW@ = Wg x , (3 = u,s). (b) the tangent planes to W@ vary continuously with 
x. (c) line elements at x are uniformly reduced in length by a factor at least C e~ Xn under 
the action of S n if tangent to or under the action of S~ n if tangent to W™, for n > 
with C, A > 0; (d) if the sign of n is changed corresponding properties hold with expansions 
replacing contractions. 

Such a system has been considered in [29] (Anosov system) and it can be shown that 
it admits a distribution p, describing the statistics of random data chosen initially with the 
Liouville distribution and that p verifies (R), being uniquely determined by the property 
stated in (R). [6] 

We note that (B),(C) imply that there are 2N — 2 positive Lyapunov exponents for S 
(a 1 ? -1 ) and that (W£) have the symmetry property W£ = iW-^. Furthermore A (and iA) 
are dense in C, although we anticipate, from the pairing property in [14] and the expected 
smoothness of the Lyapunov spectrum [22], that the Kaplan Yorke fractal dimension of A, 
[17], is ~{ 2 0(N) smaller than the dimension of C (i.e. "dimensional reduction" occurs). 
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("forcing") where u is the horizontal local velocity u = i^y (and i is the x axis unit vector) 
and coupled to a thermostat. The equations of motion are the SLLOD equations [16], [14]: 

<n = Pj/™ + mr, pj = f j - iiPyj - a Pj I 1 ) 

where j = 1, ...,iV labels the N fluid particles with mass m, Pj/m is the peculiar velocity 
of particle j, i.e. its velocity with respect to the local fluid velocity u(qj) = i7J/?; Fj is the 
interparticle force on particle j, due to a purely repulsive pair potential <p(r) where r is the 
interparticle distance (e.g. an inverse power potential). If x = (pi . . . pjvqi . . . qjv) denotes 
the phase of the system, the variable a is defined by requiring the internal energy Hq{x) = 
X/jli p|/2m + $(qi...qjv) to be a constant of motion, where $(qi . . . qjv) = J2i< 3 ^(l^ - <q_j | ) , 
and a can be easily computed [15]. The model is studied in [14], [15], [18] with periodic 
(Lees-Edwards) boundary conditions, which we replace here by simpler ones: horizontally 
periodic boundary conditions and vertically elastic wall reflections, [3]. 

For a dynamical description of this system one can suppose that the total (peculiar) 
momentum P|| and the center of mass position X|| in the shear direction vanish. If the 
observations are made in discrete time by observing subsequent particle collisions, then, with 
the constant energy, the system can be described in a phase space C of 4iV — 4 dimensions. 
The Liouville distribution can be projected on C giving a distribution fi (to which we refer 
with the same name), but fi is stationary only for 7 = 0. The evolution will be a map 
S related to the solution operator t — > S t x. of (1) and to the time t(x) elapsing between 
a timing event x £ C and the next: Sx. = S^xjx. The phase space volume variation rate, 
which is also the entropy production rate for this system [23,11] is the divergence of the 
r.h.s. of Eq.(l) and equals (2N — 2)<r, with a close to a; in fact one finds, after a brief 
computation: a = a + 7 (^1)^ ^ = a + ^fO(N~ 1 ). The average entropy production rate 
is (2N — 2)(<r) where the brackets denote a (forward) time average over an infinite time. 

The many particle statistical mechanical system in Eq.(l) exhibits the following features: 

(A) Dissipation: (a) > 0, [4]; 

(B) Time reversal invariance: the map i: (x } y } p X} p y ) — > (x } —y } —p X} p y ) is such that if 
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In his treatment of strange attractors to understand turbulence Ruelle introduced, as a 
principle (R): the time averages of observables, on motions with initial data randomly sam- 
pled with the Liouville distribution jj, 0) are described by a stationary probability distribution 
ji obtained by attributing a suitable probability density to the surface elements of the unstable 
manifolds of the points in phase space, [27]. In hyperbolic systems this leads, as a theorem, 
to identifying Jl with the SRB on the attractor, [27]. (R) was based on important previous 
results [29], [10], [26], [25], [17]. 

The unstable manifolds are quite difficult to obtain: it has been questioned whether (R) 
has any predictive value when it cannot be a priori proved. [1] It is implicit, we think, in 
Ruelle's ideas that the principle should be applied by assuming suitable properties making 
(R) valid and then retain the consequences concerning the large scale (in time and size) 
behavior. A similar situation arises in statistical mechanics where the equilibrium statistics is 
accepted on the basis of the ergodic principle (usually called hypothesis) and one derives the 
"heat theorem" (second law) and the thermodynamic properties. Analogous consequences, 
which would turn (R) into a tool for predictions, proved very hard to find in the case of 
turbulence. But recently [11] have shown that (R) does imply macroscopic consequences, 
e.g. the Einstein relation for the diffusion and conductivity coefficients in systems similar 
to ours in the linear regime where (R) can even be proven for a single particle system. 

In this letter we present consequences of a different nature, that can be derived from (R) 
and subjected to experimental tests, particularly for non small forcing, without having to 
compute rather difficult dynamical quantities like Lyapunov exponents as in [11] and without 
the need to compute, even approximately, Jl. [2] 

Although our considerations appear to be quite general, the letter is written with the 
theoretical interpretation of a numerical experiment in mind: the entropy production fluc- 
tuations in a shearing many particle fluid in a nonequilibrium stationary state far from 
equilibrium, [15]. 

We first define the model of the shearing fluid we treat. The two-dimensional shearing 
fluid is in a non equilibrium stationary state, driven by an external shear rate 7 = du x /dy 
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